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Abstract. The phase diagram for a two-dimensional self- avoiding walk model on 
the square lattice incorporating attractive short-ranged interactions between parallel 
sections of walk is derived using numerical transfer matrix techniques. The model 
displays a collapse transition. In contrast to the standard 0-point model, the transition 
is first order. The phase diagram in the full fugacity-temperature plane displays an 
additional transition line, when compared to the 0-point model, as well as a critical 
transition at finite temperature in the hamiltonian walk limit. 
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1. Introduction 

Self-avoiding walks have been widely studied as models of polymers in dilute solution|I[]. 
A self-avoiding walk on a lattice is defined as a random walk which is forbidden from 
visiting the same lattice-site more than once[|l|. In the limit of very long walks it 
models a polymer in a good solvent. In order to model the relative affinity between the 
monomers (compared with the solvent) an attractive energy is introduced between non- 
consecutively visited nearest-neighbour sitesp], ^. This is the standard 6'-point model. 
At high temperatures the polymer is happy to be in solution. As the temperature 
is lowered the polymer tends to collapse in on itself, and at a given temperature will 
precipitate from solution. For an idealised infinitely long polymer the high temperature 
and low temperature regimes are separated by a phase transition occurring at a 
temperature known as the theta temperature, Te[|], §,0. In the grand-canonical 
description the average length of the walk is controlled by a fugacity K (or chemical 
potential fi where K = exp(— At high temperatures the (average) length of the 
walk diverges smoothly as a critical fugacity Kc{T) is approached. At low temperatures 
the length jumps discontinuously at a given value of the fugacity K*(T) (corresponding 
to a first order transition). The two lines Kc{T) and K*{T) follow continuously one from 
the other, and the two behaviours are separated by a tricritical point, at a temperature 
corresponding to Tg. The transition line Kc{T) then K*{T) is identified with the 
"thermodynamic limit" of the polymer problem. Along this line, the proportion of 
visited sites is zero for T < Te and becomes non-zero continuously for T > Tq. 

While it is clear that the main experimental interest lies in three dimensions 
(although there are some experimental results also in two dimensions^), there has 
been a lot of theoretical interest in the two dimensional case[^ H, P]. This is in part 
due to the fact that the upper critical dimension for a tricritical point is three |T^. In 



dimensions lower than three (but larger than one) there is a possibility for a greater 
diversity of critical behaviours. 

In this article we investigate a two-dimensional model with a restricted set of 
interactions, compared to the theta model; an interaction is only present between 
nearest-neighbour non-consecutively visited sites belonging to straight portions of 
polymer, i.e. an interacting site cannot sit on a corner of the walk (see Figure ^). This 
choice of interactions may be motivated by the role of hydrogen bonding in the formation 
of secondary structures in proteins [0]. These bonds are formed between the CO group 



of one peptide and the NH group of a nearby (non-consecutive) peptide. This bond 
imposes a constraint on the orientation of the relevent peptides. This type of bonding 
is thought to be responsible for the formation of alpha helices and beta sheets |T2 



In our two-dimensional toy-model the formation of alpha helices is not possible. A 
representation of the formation of a beta sheet through hydrogen bonding is shown in 
Figure |l|. 

Our main interest in this paper is not to give a realistic model for protein folding, 
but rather to understand the infiuence of these apparently small modifications to the 
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Figure 1. A schematic representation of hydrogen bonds in a beta sheet (following 
Bascle et 



local attractive interactions on the critical behaviour of lattice interacting self-avoiding 
walks. To this end we present the phase diagram and elucidate the critical behaviour 
in the full fugacity — temperature plane. 

In the following section we present in detail the model as well as the transfer matrix 
calculation of the quantities of interest. In section ^ we present the phase diagram and 
discuss the different transitions present. In section ^ we concentrate more specifically 
on the Hamiltonian Walk limit of the model in which all sites are visited exactly once. 
In section ^ we finish with some concluding remarks. 

2. The model and the transfer matrix calculation 



The model studied here consists of a self-avoiding walk embedded on a square lattice. 
An interaction energy — e is assigned for each non-consecutively visited pair nearest- 
neighbour sites for which the four bonds are parallel (see figure §). A chemical potential, 
fi, and related fugacity K = exp(— are associated to each step. 
The grand canonical partition function is then 

Z=J2 K^'expime), (1) 

walks 

where N is the number of steps in the walk and Nj is the number of interactions. 
Quantities of interest may be calculated as appropriate derivatives of the partition 
function, for example, the density is given by: 
_ (AT) _ K dhgZ 
P~ Q ~ Q dK ' ^"^^ 
where Q is the number of lattice sites. The relative length fluctuations are related to 
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Figure 2. Possible configurations of nearest-neigfibour non-consecutively visited sites. 
Configurations of type (a) include an attractive interaction energy — e, configurations 
of types (b), (c) and (d) interact in tlie standard 0-point model, but here do not 
participate in hydrogen type bonding. 



the derivative of p through: 



K dp 



(3) 



In the dilute polymer regime (small enough K) the correlation function may be 
identified with the probability that two given points are joined by a walk. In this case 
the correlation length must scale in the same way as any reasonable characteristic length 
scale associated to the walk. It is usual to take the radius of gyration, which corresponds 
to the average squared distance the walk steps from the common centre of gravity. It 
follows that 



where z/ is defined as the standard correlation length exponent |TB 
(0), one finds 



(4) 

Using the relation 



Rg ~ (iv)^ (5) 

from which we may see that the walk is fractal; the mass of the walk scales linearly 
with N, and hence Rq'^ ■ The fractal (Haussdorf)|T^ dimension of the walk is then 
dff = I/'^IIH- It is usual in such problems to take equation (|^) as the definition of the 
exponent z/[^. In this case, note that it only relates to the usual critical exponent in the 
case of a critical transition in the dilute regime. In the collapsed phase(s) the polymer 
fills the lattice, the mass scales as Rq and equation gives u = 1/d. 

In this article we propose to use a transfer matrix formalism to calculate Z as a 
function of K and /3e||l5|, For convenience we shall from now on set e = 1, which 

simply corresponds to a choice of temperature scale. 
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Figure 3. Example of a polymer configuration with the corresponding column states 



The idea behind the transfer matrix formahsm is to calculate the partition function 
on a lattice which is infinite in one direction (the x direction, say) but finite in the other 
{y) direction. The thermodynamic limit is attained by increasing the lattice width. 

The standard way of considering the problem is as follows: Define the restricted 
partition function Zx{Cq, C^) as the partition function for a portion of walk between the 
origin and x. The walk has a column state Cq at the origin and in column x. One 
may then write the following recursion relation: 

Zx+i{Co,Cx+i) = ZxjCo, Cx)T{Cx, Cx+i), (6) 

where T{Cx,Cx+i) is the additional Boltzmann weight to add column x + 1 in 
configuration Cx+i next to column x in configuration Cx- This forms a (transfer) matrix. 

That this recursion should be valid for a spin system is fairly straightforward, since 
the interactions are all local. For a polymer model it is less clear that this should be 
possible, since one has to take into account non-local factors, most notably one has to 
ensure that the partition function describes only one chain, without the formation of 
"orphan" loops. This is done by appropriately defining the column states Cx- For the 
self-avoiding walk problem with no added interactions (i.e. e = 0) it is sufficient [jl5| to 
define a column configuration by the arrangement of horizontal bonds in the column 
along with information about the connectivities between the bonds, that is information 
about which pairs of horizontal bonds are connected by polymer loops to the left (taking 
X as increasing towards the right). See figure |] for an example. By successive application 
of (H), one finds 

Zx{Co,C.) = {Co\r^\Cx). (7) 

The situation becomes a little more complicated if interactions are to be taken 
into account. It is necessary to know the column states over three columns in order to 
incorporate the horizontal bonds. The appropriate modification to the above equations 
is 



16, 18 



Zx+i{Co,Ci;Cx,Cx+i) — ^ Zx{Cq,Ci;C'x-iC'x)T[C'x-i,C'x',Cx,Cx+i),{8) 
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and 

ZxiCo,Ci;Cx-i,Cx) = {Co,Ci\T^ ^\Cx-i,Cx)- (9) 

The primed configurations are introduced to retain the matrix notation, the elements 
of the matrix are now labelled by the two column states for the input and two for the 
output (primed), with the restriction that only matrix elements with Cx = C'x niay be 
non-zero. 

It is required that the partition function sum over walks of all possible lengths. To 
do this we define a restricted partition function Zx for all the walks having a caliper 
extension x, i.e. walks which are entirely confined to a strip of length x, and extend 
over the entire length of the strip. This is simply Zx summed over all configurations 
Co, Ci, Cx-i and Cx compatible with the starting and ending configurations of the walk. 
The full partition function is then given by 

oo 

z=Y.Zx+^, (10) 

x=0 

oo 

= E'E^^ (11) 

b.c. x=0 

where the primed sum is over the boundary configurations, as described above. 

Let us denote by {Aj} the eigenvalues of T, numbered in order of their moduli. The 
eigenvalue of largest modulus, Ai , is always positive. 

For small enough i^', Ai < 1. In this case the partition function may be written 

a 

The ttj are essentially constants dependent on the boundary conditions, which gives for 
the density: 

If the derivative of A, is finite, which is the case at least for small enough (3, the density 
is zero as long as Ai < 1 and may only change from its zero value when the average 
length of the walk diverges, i.e. Ai > ifl^, |18|. The point Ai = 1 is identified with the 
critical transition, and is often used in polymer science to define the "thermodynamic 
limit" of "very long" polymers[|l|. Another way of reaching to the same conclusion is 
to consider two points as being correlated if they are joined by a walk. One may then 
identify Zx with a two point correlation function, and hence identify the correlation 
length, ^0: 



log ( AT j 

The correlation length diverges as Ai — ^ 1, consistent with the identification of a critical 
transition. 

The above argument needs to be treated with caution. As long as the lattice width 
is finite, the system is equivalent to a one-dimensional system with a finite, if large. 
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number of states per site (each site corresponding to a column in our model). It is well 
known that a one dimensional system with a finite number of states per site cannot have 
an equilibrium phase transition [1^ . 

The reason for this apparent discrepancy is that the definition of ^ ( p!^ is not valid 
in the large K phase (dense polymer phase). In a spin model the partition function 
would have been given by in the limit x —>■ oo. The transfer matrix would contain 
all the possible configurations, and a standard analysis shows that [l20| : 



(15) 

where Aq is the largest eigenvalue, and Ai is the second largest (in modulus). To ensure 
the correct construction of a single walk all configurations consisting of an empty column 
have been excluded from the transfer matrix, but replaced by the sum over all possible 
polymer extensions x. The consequence of this choice is that the largest eigenvalue in 
the zero-density phase has been eliminated. Its value, as can be seen by inspection of 
the expression for ^, would be Aq = 1. 

This prescription only holds as long as Ai < 1. If Ai > 1 then the partition function, 
density and ^ are given by: 

Z ~ lim Af , (16) 
P = LX[dK ^''^ 

e = , (18) 

where L is the lattice width. 

For convenience we impose the condition that the walk of extension x' in the x 
direction originates at x = and terminates either at a; = or x = x'. As is 
usual in equilibrium phase transitions, the precise choice of boundary condition does 
not change the critical behaviour. Imposing this restriction it is easy to see that 
column configurations with even numbers of horizontal bonds may only follow column 
configurations with even numbers of horizontal bonds, and likewise for configurations 
with an odd number of horizontal bonds. This block diagonalisation is convenient since 
it transpires that for much of the phase diagram Ai and A2 arise in different blocks. 
It is relatively easy to calculate the largest eigenvalues of a matrix using the power 
met ho dpi . 



3. The K — 13 phase diagram 

The critical lines may be identified using phenomenological renormalisation . At a 
critical point the correlation lengths for two strip widths, measured as a fraction of the 
strip width, must be the same. This refiects the scale invariance of a critical system. In 
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practice a finite width system is always off-critical, however 



(19) 



L L' 

gives a finite-size estimate of the critical fugacity, K* , which will tend to the true critical 
fugacity as the strip widths tend to infinity. Assuming that the system is sufficiently 
close to the true critical point, the correlation length behaves (to leading order) as 

iL = A\K*-K,\~''. (20) 

Since it is the finite width of the system which prevents the correlation length from 
diverging, at K* we have L oc ^l- If these two scaling laws are admitted, then a finite- 
size estimate of the correlation length exponent may be calculated from the equation: 

— - 7";lr pi) 

"i.i' log (i) 

There exists at least one phase transition. In the case f3 = this corresponds to 
the standard Self- Avoiding Walk transition. As j3 is increased, this point extends into a 
line of critical points. In analogy to the standard 6'-point model, we would expect this 
transition to change to first order at some given value of /3 = Ph- As discussed in the 
previous section, starting from the low-K phase, the correlation function is given by ([T^). 
In the high-K phases the appropriate form of is given by (|15]). The phase diagram 
estimates found using phenomenological renormalisation, using the appropriate form of 
^, are shown in figure ^. As is usual in this sort of problem, there are strong parity 
effects, notably for the odd lattice widths Ai is in the odd sector of the transfer matrix, 
while for the even lattice widths Ai is taken from the even sector. The next largest 
eigenvalue, A2, is in the other sector for the low-i^' phase and for the high temperature 
high-i^ phase, while it is in the same sector for the high-ii" low temperature phase. 
This observation was not used to estimate the transition lines, which were however 
found assuming that Ai and A2 came from different sectors of the transfer matrix. Not 
all the transition lines are found for even lattice widths. 

Let us now characterise the different phases and transitions. For low K the average 
length of polymer is finite, the density of the polymer on the lattice must therefore be 
zero. The first transition as K is increased is to one of two dense phases, depending on 
whether P < Ph or (3 > jSn- It is convenient to use the density as an order parameter. 
For practical reasons we choose boundary conditions in which the polymer is made to 
extend the length of the lattice, between x = and x 00. This obliges the polymer 
to be infinite in length even in the zero-density phase. A choice of boundary conditions 
is not expected to affect the critical behaviour of an equilibrium system. The partition 
function becomes Z = lima;^oo Af for both the high and low density phases, and p may 
be calculated using ( p!7D in all phases. The density will now be finite for K < Kc as long 
as the lattice width is finite, but must tend to zero as L ^ 00. This is indeed what is 
observed from the density plots shown in figures |^ and ^. 

For P < Ph the behaviour of p and the length fiucuations indicate the existence 
of a line of second order phase transitions, as expected. The line is in the self-avoiding 
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Figure 4. Phase diagram estimates found using phenomenological renormalisation 
group for odd lattice widths (top) and even lattice widths (bottom) 

walk universality class, having = 3/4 as may be seen from the finite size estimates for 
u shown in figure |^. The transition is to a finite density phase. 

For (3 > (3h the density plot, given in figure ^, indicates that the transition is first 
order to a dense phase. That this line corresponds to the compact polymer state is 
confirmed by the observed exponent v = 0.5 = 1/d, see figure |^. A major difference 
with the ^-point model, seen using the phenomenological renormalisation scheme is the 
existence of two different high density phases. The density plots in figure ^ for /5 > Ph 
suggest that the walk changes discontinuously from p = to p = 1, and that the 
density of interactions saturates instantly. In other words the lattice changes abruptly 
from being essentially empty to essentially full, with a maximum density of interactions. 
Trivially, a lattice which is empty (a finite walk length on an infinite lattice) has zero 
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Figure 5. Plots of density, p (top), and density fluctuations (bottom) for the inverse 
temperature /3 = 0.7 
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Figure 6. Plot of density p for /3 = 1.1 



1.5 - 



0.75 



0.5 - 



0.25 - 



sizes 3 and 5 

sizes 5 and 7 

sizes 7 and 9 




0.2 0.4 0.6 „ O.f 



I 1.2 



Figure 7. Estimates for tlie exponent v for tlie low- if transition line 



energy and zero entropy per site; the free energy per site is zero. If the lattice is 
completely full with the interactions density saturated, the walk may only be in one of 
two configurations; either all the bonds are vertical or all the bonds are horizontal. The 
entropy per site in this phase is therefore also zero. Since, in this state, there is one bond 
and one interaction per site, the energy per site is e = /i — e. The transition between 
the two phases would then occur when the free energy per site, / = e — Ts = e = 0, in 
other words when /i = e. Remembering the choice e = 1, the transition line would be 
K = exp{—j3). This line has been plotted for comparison on the phase diagrams, shown 
in figure The agreement with the observed first-order transition line is excellent. 




Figure 8. Plots of the density p (top), the specific heat C (middle) and the density 
of interactions pi (bottom) for iiT = 1.5 
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In order to differentiate between the two high density phases, it is convenient to use 
the interaction density, pj, as an order parameter. This order parameter, along with its 
ffuctuations, is plotted as a function of /3 for = 1.5 in figure |[ The order parameter 
differentiates between the two phases since it is identically one in the high-/? phase, 
and less than one in the low-/? phase. The variation of the order parameter appears to 
be fairly smooth for all lattice sizes, and leads one to conjecture that this high density 
transition is critical. Also plotted in figure || is the specific heat, calculated as the second 
derivative of the free energy with respect to the temperature, keeping /i and e constant. 

Various estimates of the critical (3 ioi K = 1.5, calculated using the 
phenomenological renormalisation group and by locating the peaks in the specific heat, 
are given in table |I|. The various results are coherent with an estimate l3c{K = 1.5) = 
.55±0.01. At a critical transition, the specific heat diverges as T —s> Tc with an exponent 
a, 

C~|T-Te|". (22) 
When combined with relation pO| ) and oc L, this leads to the finite-size scaling result 

C^ax~A + 5L"/^ (23) 

A three point fit for odd widths gives ajv = 0.67 for L = 3,5, 7 and ajv = 0.55 for 
L = 5, 7, 9. A similar three point fit for the even sizes, L = 4, 6, 8 gives a value of 
ajv = 0.156. If sufficiently large lattice widths are considered then A can be dropped. 
When this is not the case, then the result is an effective value of the exponent, which 
should converge to the correct value as the lattice width is increased. Unfortunately the 
number of two point approximations we could construct does not allow for meaningful 
extrapolation. All we could determine is that 0.156 < a/v < 0.55 gives a reasonable 
limit on the possible values of a/ v. Using the scaling relation a = 2 — du, this gives u 
as 0.78 < V < 0.93. Note that this range is coherent with the values of v found from the 
phenomenological renormalisation group calculation, shown in table |I[ Additionally, 
the upper boundary, calculated from the even lattice sizes, is in good agreement with 
the even lattice size estimates of z/ in table I. 

It now remains to determine the nature of the transition at {Kh,Ph), where the 
three transition lines join. For the standard 6 point model, this point is tricritical, and 
corresponds to the confluent of two transition lines, not three. 

The problem we are faced with is that the phenomenological renormalisation group 
applied to even lattice widths does not give the first-order transition at all. When 
applied to the odd lattice widths, the method is not able to follow the self-avoiding walk 
transition line all the way to Ph, jumping to the upper transition line just short of Ph- 
This undershoot rapidly becomes smaller as the lattice width is increased. Remembering 
that the length of the walk first diverges when Ai ^ 1 and that this condition coincides 
with the low-K phase transitions in the thermodynamic limit, we can use the condition 
Xi{K,l3) = 1 as a criterion for defining finite width estimates to the transition lines. 
The phase diagram estimates calculated using this criterion are shown in figure ^ This 
method is incapable of giving the high-i^' transitions, and converges more slowly than 
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Table 1. Estimates for (a) /3c and v for K = 1.5 limit using phcnomcnofogical 
renormalisation group with strips of width L and L + 2 and (b) Pc as estimated from 
the peak of the specific heat and the height of the peak 

(b) 



L 




1/ 




C'max 


3 


0.579774 


0.885482 


0.572447 


0.920518 


4 


0.546148 


0.920203 


0.500045 


0.645603 


5 


0.560960 


0.841050 


0.561950 


1.344506 


6 


0.550013 


0.920194 


0.510141 


1.131212 


7 


0.554531 


0.864167 


0.556259 


1.706669 


8 






0.524612 


1.494860 


9 






0.553865 


2.024505 



the phenomenological renormalisation group, but has the advantage that the transition 
hues are fairly smooth and continuous. The density and its fluctuations, calculated 
along this line, are presented in figure |10|. As the width of the lattice is increased, 
the jump in the density becomes sharper. The height of the fluctuations (shown with 
a log-linear scalin in figure |rT]) increases extreemly rapidly with lattice width, stronly 
suggesting a delta peak in the thermodynamic limit. Whilst the range of lattice widths 
prevent us from being categorical, the results suggest strongly that the transition is first 
order along the transition line. This seems to confirmed in the phase diagram. In both 
figures § and | the transition line has a discontinuous slope at Ph- 



The positions of the peaks of density fluctuations are shown in figure [IT] along with 
the positions of the cusps where the upper and lower critical lines join for even lattice 
widths. We estimate Ph = 1-00 ± 0.02, with = exp(— /?/^). It is curious that Ph 
should be so close to one, though we have as yet not been able to determine why this 
should be the case. 

4. The Hamiltonian Walk limit (K oo) 

In the limit K —* oo the walk will fill the lattice maximally. Each site will be visited 
once and only once. In this limit the walk becomes what is known as an Hamiltonian 
Walk[p^. In this limit the effect of the interaction is to energetically penalise corners 
in the walk. This is similar to the Flory model [^] except that there are configurations 
which introduce one or two corners into the walk with no difference in energy, see 
figure 0. 

There is a low-temperature corner-free phase. As the temperature is raised there is 
a transition in which the density of corners becomes non-zero. To distinguish between 
the high and low temperature phases it is possible to use either the density of corners 
the energy density. This latter is shown, as a function of /? in figure |l^. 

We estimate the critical temperature again using the phenomenological 
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Figure 9. Phase diagram estimates calculated using the condition that the largest 
eignvalue Ai = 1 odd lattice widths (top) and even lattice widths (bottom) 



renormalisation group, the results of which are shown in table || and extrapolated using 
a three point fit to 

13^ = (3^ + AL^, (24) 

where I3c is the finite size estimate and is the extrapolated value. A and B are 
constants to be determined. 

Only the odd lattice widths presented solutions to the phenomenological 
renormalisation group equation (p!9D, given in table |^. The positions and values of 



the maxima of the energy fiuctuations per site, shown in figure |T3|, are also reported in 
table 0. The fiuctuations appear to saturate. We assume that, even if the fiuctuations do 
not diverge, the transition point is estimated by the value of (3 for which the fiuctuations 



Two dimensional self- avoiding walk with hydrogen-like bonding 



16 



A 
V 



A 

•z. 

V 

I 

A 
< 

V 




Figure 10. Density (top) and density fluctuations (bottom) along the transition line 
as estimated using the condition that the largest eigenvalue Ai = 1 



are maximal. This agrees well with the phenomenological renormalisation group 
estimate, and fits with a continuation of the finite K phase diagrams (figure ^ to high 
values of K. Phenomenological renormalisation and the peaks of the energy fluctuations 
give estimates of the critical temperature which concur, giving f3c = 0.476 ± 0.001. 

It is also possible to estimate the exponant r], using the conformal invariance result 
for periodic boundary conditions[^]: 



L 



We find r]^ = 0.096 ± 0.002, see table |. 

It is anticipated that the entire high temperature phase is critical, in analogy with 



other models of this type, such as the model F|^, the Flory model |^ and the loop 



Two dimensional self- avoiding walk with hydrogen-like bonding 



17 



0.9 - 



0.85 - 




Figure 11. Finite size estimates of Ph, the value of /? at the collapse transition. The 
upper points, calculated for even lattice widths, come from the estimated position at 
which the two critical lines merge. The lower points, calculated for the odd lattice 
sizes, come from the peaks of the fluctuations in N. The error bars given correspond 
to the interval between the points calculated. These error bars could be reduced at an 
additional cost in computer time. 



Table 2. Estimates for (a) /3c, ly and r] in the K ^ oo limit using phenomenological 

rcnormalisation group with strips of width L and L + 2 and (b) (3^ , the critical value 
of (3 as estimated from the peak of the interaction fluctuations, the height of the peak 
(./max) and the value of t] calculated the extrapolated value of f3c- Even and odd lattice 

widths arc^ lakt^u ajjart 

(5 (b) 



L 




1/ 


V 


f3^ 






3 


0.450472 


1.109084 


0.082484 


0.435221 


1.714121 


0.087226 


4 








0.375742 


1.875099 


0.035357 


5 


0.461359 


1.060628 


0.085672 


0.465233 


2.124762 


0.090245 


6 








0.437435 


2.260666 


0.047824 


7 


0.465970 


1.105916 


0.087611 


0.471160 


2.354970 


0.092053 


8 








0.454834 


2.491249 


0.056830 


9 








0.473274 


2.492086 


0.093240 


10 








0.461866 


2.636003 


0.064064 


L — > 00 


0.477 






0.476 ±0.001 


2.867 


0.096 ± 0.002 


(odd) 














L — > 00 
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Figure 12. Energy density as a function of (3 in the K oo limit for odd lattice 
widths (top) and even lattice widths (bottom) 



gas model on the brickwork lattice |2^. The estimates of u are not at all clear. They 
seem close to u = 1, but in other respects the transition is similar to that seen in the 
model F and the Flory model, which have transitions of infinite order with u infinite. 



5. Conclusions 



In this article we have presented a two-dimensional model for homopolymer collapse 
under the influence of "hydrogen-bonding" like nearest-neighbour interactions. The 
phase diagram is quite different from the equivalent 6 point phase diagram; the collapse 
transition is now first order. This is also observed in the oriented interacting self- 
avoiding walk model, where, for appropriate values of the model parameters, a first 
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Figure 13. The energy fluctuations in the K ^ (x hmit for odd lattice widths (top) 
and even lattice widths (bottom) 



order collapse occurs to a spiral configuration! 28 1. These two models have in common 



that the collapsed phase is anisotropic, in our case one of the two lattice directions is 
selected, and in the oriented walk case, for strong enough parallel interaction strength, 
one of two chiralities is selected. In the standard 9 model the collapsed phase is isotropic. 

Here we also observe an additional transition in the dense walk regime. This 
transition appears to be second order, though for the moment we have not been able 
to determine to any accuracy any of the critical exponents. Similar critical lines 
are observed in the loop models with corner interactions, with or without a collapse 
transition [p^, p^ . 

There also exists non-trivial critical behaviour in the K oo limit, where there is 
a finite temperature phase transition between a critical high temperature phase and a 
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frozen low temperature phase. 

Another well studied model with a qualitatively similar phase diagram is the 



partially directed interacting polymer model[30, M. There are major differences 



however: The line K = exp(— /?) corresponds to a phase transition line for all values 
of K, deliminating a similar fully collapsed phase. The equivalent phase to the finite-p 
phase has a density p = due to the directed nature of the model. The equivelent to 
the theta point transition is second order |Q (here it appears to be first order). 
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